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A two-leg ladder with n-component fermionic fields in the chains has been considered using an 
analytic renormalization group method. The fixed points and possible phases have been determined 
for generic filling as well as for a half-filled system and for the case when one of the subbands is 
half filled. A weak-coupling Luttinger- liquid phase and several strong-coupling gapped phases have 
been found. In the Luttinger liquid phase, for the most general spin dependence of the couplings, 
all 2n modes have different velocities if the interband scattering processes are scaled out, while n 
doubly degenerate modes appear if the interband scattering processes remain finite. The role of 
backward-scattering, charge-transfer and umklapp processes has been analysed using their bosonic 
form and the possible phases are characterized by the number of gapless modes. As a special case 
the SU(n) symmetric Hubbard ladder has been investigated numerically. It was found that this 
model does not scale to the Luttinger liquid fixed point. Even for generic filling gaps open up in 
the spectrum of the spin or charge modes, and the system is always insulator in the presence of 
umklapp processes. 



PACS numbers; 71.30.-|-h, 71.10.Fd 



I. INTRODUCTION 

Recently, systems exhibiting non-Fermi-liquid behav- 
ior have been intensively studied, especially in connec- 
tion with the unusual normal-state properties of high- 
temperature superconductors. Low-dimensional models 
are of special interest in this respect due to their possible 
Luttinger-liquid^ behavior. The simplest model that has 
a non-Fermi-liquid ground state is the exactly solvable 
one-dimensional (ID) Hubbard chaini^ Unless the band 
is half filled, this system is a prototype of Luttinger liq- 
uids. 

In most strongly correlated systems, the degenerate d 
or / bands of the transition-metal or rare-earth ions play 
an important role and the orbital degrees of freedom, too, 
have to be taken into account. In an obvious although 
highly simplified generalization of the Hubbard model, it 
is assumed that the spin and orbital degrees of freedom 
can be treated as a unified degree of freedom with n pos- 
sible values. The SU(n) symmetric generalization of the 
SU(2) symmetric ID Hubbard model has been studied by 
several authors^di^ to check its possible Luttinger-liquid 
behavior and to study the metal-insulator transition at 
half filling or 1/n filling. It has been found that the be- 
havior is qualitatively different for n > 2 and n = 2. 

The physics is more complicated when — as a first step 
towards two-dimensional systems — two Hubbard chains 
or equivalently two Luttinger liquids are coupled by one- 
particle hopping into a ladder. The first question that has 
been addressed is the relevance or irrelevance of the in- 
terchain hopping. It has been shown^'^'^ that for generic 
filling the ID Luttinger liquid is unstable with respect to 
arbitrarily weak transverse single-particle hopping. Ei- 
ther the system becomes Fermi liquid or new types of 
states may occur. The situation may be different at a 
special filling where umklapp processes may suppress the 
single-electron interchain hopping^ 



A more difficult problem is to decide what these states 
may be. To answer this question the Hubbard ladder — or 
more general models with relevant intrachain backward 
scattering processes, too — have been studied by differ- 
ent analytic and numerical methods like renormalization 
group in bosoniSiiSiiiiiSiiii^ or fermioni^*i& representa- 
tion, Monte Carlo simulationsiii^ or using the density- 
matrix renormalization group ^SiSSiSiiSS 

As a possible classification of the new states Balents 
and Fisher— proposed to characterize the phases by the 
number of gapless charge and spin modes. Since a two- 
leg ladder with two-component fermions on both legs has 
two charge and two spin modes, in principle nine phases 
can be distinguished in this way. They have shown that 
seven of them may appear in the weak-coupling regime of 
the Hubbard ladder. The others can be realized in more 
general models only, which are defined by more coupling 
constants. 

A more complete characterization of the phases can 
be obtained by investigating the appropriate correlation 
functions. When only forward scattering processes are 
allowed for on the chains, for generic filling, the domi- 
nant singularity appears in the spin-density and charge- 
density responsesfii while if backward-scattering pro- 
cesses are also taken into account the two-chain model 
has predominant c?-type pairing fluctuations even for 
purely repulsive interactions. ^"'^^ Similar result, the pos- 
sibility of density wave or superconductivity has been 
found in a variational calculation.— 

Similar situation occurs for the half-filled Hubbard lad- 
der. This model has been shown^- to scale to the S0(8) 
Gross-Neveu model which has fully gapped low-lying ex- 
citation spectrum. For repulsive interaction the system 
is a Mott-insulator spin liquid with d-wave pairing cor- 
relations. It is interesting to note that if both repul- 
sive and attractive interactions are allowed for, the phase 
diagram contains four phases related by SO (5) symme- 
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try, which can be interpreted^ as unifying magnetism 
and superconductivity. Similar approximate SO (5) sym- 
metry has been foundS^ at low energies in the repul- 
sive strong-coupling regime of two coupled Tomonaga- 
Luttinger chains. 

The competition of different density-wave or supercon- 
ducting orders has been studiedi^^ in a more general 
model, where besides the one-particle transverse hopping 
the Coulomb interaction and Heisenberg exchange be- 
tween sites on the same rungs and along the diagonals of 
the elementrary plaquette as well have been taken into 
account. 

Although the problem of two coupled two-component 
fermionic chains has been extensively studied, and some 
results on the two-dimensional SU(7i) Hubbard model are 
knownjSiSl to the best of our knowledge the properties of 
two coupled SU(n) symmetric fermionic chains have not 
been discussed. For this reason, in this paper, we will 
present the results obtained for a ladder built up from 
two n-component fermion systems. In this respect this 
calculaton is an extension of our earlier worki^ On the 
other hand, since a two-leg ladder can be treated as a 
two-band model, and the method applied is the multi- 
plicative renormalization group combined with bosoniza- 
tion treatment, this paper can also be considered as a 
generalization of Refs. 16 and to n-valued spin. 

The model is introduced and the types of interac- 
tions are discussed in Sec. ^ The next section contains 
the scaling equations. The behavior of the system in 
the weak-coupling Tomonaga-Luttinger fixed point and 
the properties of the gapless Tomonaga-Luttinger phase 
of the two-band model are analyzed in Sec. IV. Here 
the calculations are extended to fully anisotropic spin- 
dependent couplings. The role of the backward, charge- 
transfer and umklapp scattering processes is discussed in 
Sec. V, and the possible gapped phases characterized by 
the number of gapped modes are presented. The numeri- 
cal results obtained for the Hubbard ladder and for some 
more general models are presented in Sec. VI. Finally 
Sec. VII contains a summary of the main results. 



with j = 1,2, are coupled into a ladder by a rung hopping 
term 

N n 

^ E E («M,.«2,^,. + 4,.,.ai,.,.) • (2) 

4=1 CT=1 

The local U(n) symmetry of the interaction term is re- 
duced to global U(n) symmetry by the kinetic energy 
term. An SU(n) symmetric model is obtained if the over- 
all U(l) phase factor is removed. Although in real space, 
the Hamiltonian has a particularly simple form, in what 
follows it will be treated in momentum representation, 
where it is easier to separate the relevant, irrelevant or 
marginal components. 

A. Kinetic energy 

The terms describing the intrachain hopping are diag- 
onal in momentum representation. When the interchain 
hopping is taken into account, the full kinetic energy is 
diagonalized by the symmetric and antisymmetric com- 
binations of the operators belonging to the two chains: 

bk,<j — —^{ai,k,<7 — a2,k,a)- 

As it is shown in Fig. ^ the spectrum consists of two 
bands labelled by a and b which are obtained by shifting 
the dispersion curve of the free chain by ±tj_. 
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II. THE MODEL 

The simplest model we consider is a two-leg SU(n) 
Hubbard ladder in which two identical SU(n) Hubbard 
chains described by the Hamiltonians 

N n 



^^=^EE(«k.«..m.. + H.c.) 



u ^ 



(1) 



^EE 

i—l 17,(7' — 1 



FIG. 1: The spectrum of two coupled chains: the spectrum 
before (after) hybridization is shown by dotted (dashed) lines. 
Solid lines indicate the linearized spectrum. 

When the number of electrons is small, the upper band 
is completely empty, and the behavior is identical to that 
of a one-band Hubbard model. The behavior is similar 
in the case when the system is close to saturation, and 
the lower band is completely filled. In what follows it 
will be assumed that the band filling is between these 
two extremes, both bands are partially filled, moreover, 
the Fermi momenta are not too close to the zone center 
or zone boundary. Under these circumstances one can 
assume — as is done usually — that the relevant electronic 
states are the ones near the Fermi points and their dis- 
persion relation can be approximated by straight lines 
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with slopes ±hva and ±hvb around the Fermi points ±fca 
and ±fcb, respectively. If the lattice constant is taken to 
be unity 



fi'Va(b) — 2t sin k, 



■a(b) 



(4) 



and the value of the Fermi momenta is determined by 
the overall filling of the band and the splitting of the 
subbands, 2t±. It follows from the symmetry of the 
tight-binding dispersion relation of the free chain that — 
as shown in Fig.|21— for a half- filled system Va = Vb, while 
more generally Va 7^ ffc. 



2t 
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FIG. 2: The linearized spectrum for a half-filled ladder. 

Measuring the momenta from the respective Fermi 
points the kinetic energy takes the form 



^kii 



J2 [^«afc(ai„ 



-ka+k,<y"'-ka.+k,(7) 



^^&^(^fcb + fc,cr^fci,+fc,cr ^-kf, + kM^-kf, + kM ) 



(5) 



The operators al^^+fc^,, («±fc„+fc,^) and 6^^^+^^^ 
i^±ki+k a) create (annihilate) an electron with momen- 
tum k and spin a around the Fermi points ±ka and ±fcb, 
respectively. Thus, the linearized spectrum consists of 
four branches, two corresponding to right-moving parti- 
cles and two to left movers. In what follows the notation 



'^-.k.cr — ^-ka + k,a- ' 
^+,k,<7 = ^kb+k,a ' ^-,k,<7 = b-kb+k,cr 



(6) 



will be used. 

For simplicity a sharp cut-off is assumed in the allowed 
momenta, |fc| < fco, leading to bandwidths 2Ea{b) = 
2hVa(b)ko, which in general are different. It is also as- 
sumed that the bandwidth is comparable or less than 
the splitting of the bands, in other words the difference 
ka ~ kb is comparable or larger than the cut-off fcg . Due 
to this choice, from the results presented in this paper, it 
is not possible to recover the limit of uncoupled chains. 



B. Interactions 

When the on-site Hubbard interaction is written in 
terms of the operators of right and left-moving fermions 



of type a and 6, the terms appearing have the generic 
form 



,k4,cr ' 



(7) 



where the operators , , 7 and S are creation and anni- 
hilation operators of particles of type a or b, and Ai = ±. 
While in the Hubbard ladder the strength of all these 
processes is related to the single on-site Coulomb repul- 
sion U , in what follows a more general Hamiltonian will 
be used allowing for different couplings for the different 
scattering processes. 

The coupling constant of the process given in 10 will 
be assumed to be momentum independent and will be 
denoted by 



It is clear that the terms with couplings 

and 



af3yS 
yAi,A2,A3,A4,crcr' 



(8) 



(9) 



describe the same process. Similarly, the terms with cou- 
plings 



•y \i , A2 , A3 



, A4 jCTcr' 



and 



57/3 ct 

^A4,A3,A2 



. Ai .aa' 



(10) 



are Hermitian conjugates to each other. Assuming that 
the couplings are real, the same coupling constant is used 
for them. If the system has inversion symmetry, the cou- 
plings 



^Ai , A2 , A3 ,A4 ,0 



and 



Q/3715 

^— Ai, — A2, — A3, — A4 ^aa' ' 



or when combining with (jSJ the couplings 



a/375 



and 
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-Ai, — A4, — A3, (t'ct 



(11) 



(12) 



should be taken to be identical. Furthermore, since we 
have taken the symmetric and antisymmetric combina- 
tions of the operators belonging to the chains, only such 
scattering processes survive in which all four particles are 
in band a or in band 5, or two of them are of type a and 
two of them are of type b. 

A scattering process is intraband if all participating 
particles belong to the same subband (a or 6), a = /3 = 
7 = (5. The processes in which either a (3 or a 6 are 
interband. The terms with coupling of type (jr""''^ describe 
processes in which charge is transferred from one band to 
the other. If a =/= cr' in an interband process with coupling 
gabab^ the Scattering is accompanied by a spin flip in the 
bands. This could be considered as the analogue of the 
Kondo coupling. 

It can be assumed that in a non-polarized system the 
coupling constants depend on the relative orientation of 
the spins only. As a further constraint we will distinguish 
two cases only, whether a = a' or a ^ a' . The index | 
or _L will be used for them. 

Momentum conservation gives two further restrictions: 



fci 



(13) 
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and 

Aifca + Aafc^ = Agfe^, + + G, (14) 

where G is a reciprocal lattice vector. G — corre- 
sponds to normal processes, while G is finite for um- 
klapp processes. Since it was assumed that ka and kf, are 
sufficiently different, momentum conservation (including 
umklapp processes) can be satisfied only if: 

1. All four particles participating in a scattering pro- 
cess are right movers or all of them are left movers, 

Ai = As - A3 = A4 . (15) 

Conventionally these processes are labelled by an 
index 4. They will be neglected since — as opposed 
to the other processes — their contribution is not 
logarithmically divergent in perturbation theory. 
Normally they lead to a Fermi velocity renormal- 
ization only. 

2. Two right (left) movers are scattered into two left 
(right) movers, 

Ai = A2 = -A3 = -A4 . (16) 

These are the umklapp processes and they will be 
labelled by an index 3. 

3. A right and a left mover scatter on each other and 
again a right and left mover are created. Depending 
on whether the momentum transfer between the 
particles of the same spin is small 

Ai = -A2 = -As = A4 (17) 

or large, 

Ai = -A2 = A3 = -A4 (18) 
the processes are labelled by indices 2 or 1. 

The momentum transfer is small, the process is forward 
scattering, if either fc^ — fc^ = (and consequently kj — 
kp = 0) or ka — kg — ±(fca — fcf,). As shown in Fig. |31six 
such processes are possible. 

Besides the intraband processes, either in band a or 
in band 6, there are four possible intcrband processes. 
Two particles from band h are scattered into band a, or 
vice versa, in channels a^a^bb and b^b^aa, and a small 
momentum of order ±{ka ~ kh) is transferred between 
right-moving and left-moving particles, respectively. The 
transferred momentum can be arbitrary small in the pro- 
cesses a^b^ba and b^a^ab. An a^b^ab process with small 
momentum transfer is possible only when all particles are 
right or left movers, and these processes — as mentioned 
above — will be neglected. 

Taking into account that the processes a^a^bb and 
aa are Hermitian conjugates and the process b^a^ab 
can be obtained from an a^b^ba process by inversion, the 
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FIG. 3: Forward ((72) scattering processes in a ladder system. 

forward-scattering processes are characterized by 8 cou- 
pling constants: 
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(19) 



For six of these eight forward-scattering processes, the 
particles remain in the same branch of the spectrum. As 
we will see, they play a special role. Due to the conser- 
vation of the number of particles and of the spin in each 
branch, if they act alone, Luttinger-liquid-like behavior 
is found. 

In the large-momentum-transfer processes, Ai = — A4 
and A2 — —A3, and consequently k^ — ks = ±2fca, ±2fcf, or 
±(fca -I- kb), one has to distinguish normal (labelled by 1) 
and umklapp (labelled by 3) processes. Considering first 
the normal backward-scattering processes, two of them 
{d^a^aa and b^b^bb) correspond to intraband processes, 
and four of them [a^a^bb, 6^6^aa, a^b^ab, and 6^a^6a) are 
interband processes. The processes a^b^ba and b^a^ab are 
forbidden by momentum conservation for finite hopping 
amplitude. The six allowed normal backward-scattering 
processes are shown in Fig. 0] 

Since the processes a^a^bb and b^b^aa are Hermi- 
tian conjugates to each other and the processes a^b^ab 
and b^a^ba are related by inversion, these processes — 
analogously to the forward-scattering ones — could be 
characterized by 8 coupling constants. The backward 
and forward-scattering terms are, however, not indepen- 
dent when the particles have identical spins. Looking at 
the terms with couplings 5°^"^ and g^^°' given by 

}_ ST^ ( abab t u 
T V9l|| "'+,k+q,a°-,k'-q,a"-+.k'.a°-,k,c! 

k,k',q (20) 

+ 92t^U+qA.k'^q,.b^.,,^„a^,,J, 
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FIG. 4: Backward (gi) scattering processes in a ladder sys- 
tem. 
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FIG. 5: Umklapp ((73) processes in a ladder system. 



these are the processes in the last rows of Figs. |31 and 
01 it is easily seen that — due to the indistinguishability 
of the fermions — they describe exactly the same process, 
if the model is defined by a bandwidth cut-off, as is the 
case in our model. There is only a sign difference in 
their contribution. Similar arguments are valid for the 
couplings 35*11 and and also for gj* and g2|| • For this 
reason we use the combinations: 
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terband umklapp processes with different coupling con- 
stants could be distinguished. 

In order to get simpler formulas in what follows the 
dimensionless couplings 



and the quantity 



9. 



^-Kh{Va + Vf3 + + Vs) 



1 



^VaVb 



(23) 



(24) 



will be used. 7 characterizes the relative slopes of the 
two bands. 7 = 1 if and only if Va = Vb, while 7 > 1 
otherwise, if both Va and Vb are positive. 



and only the processes with couplings g'^_^ , g\_^ , , and 
gl*^ are true backward-scattering processes. 

When either band a, band b or the whole spectrum is 
half- filled, it is possible that Xika + X2kfs = X^k^ + X^kg + 
G with a vector of the reciprocal lattice. In the first case 
intraband umklapp processes are allowed for, while in the 
second case there are four kinds of interband umklapp 
processes as shown in Fig. [S] but as mentioned above, 
two of them appear with the same coupling constant. 

Due to the requirement of antisymmetry of the 
products of fermion operators, no intraband umklapp 
processes are allowed between particles of the same 
spin. This is also true for the charge-transfer umk- 
lapp processes of type a\ 



However, in- 
b_ „ and 



4. 

terband umklapp processes of type a\_ ^b^ ^ 

a!!^_g.6^_ £ra_ £r are not excluded by symmetry. Since 
they correspond to the same type of process the notation 



ab abab 

i/sll — il3\\ 



abba 
93\\ 



(22) 



will be used. Thus altogether two intraband and four in- 



III. RENORMALIZATION-GROUP 
TREATMENT 

The perturbative corrections to the vertices and self en- 
ergy of the ID interacting fermion system defined by the 
Hamiltonian in and ||2Jl are logarithmically singular. 
A convenient procedure to sum up the leading or next-to- 
leading logarithmic corrections to any order of the cou- 
plings is provided by the multiplicative renormalization- 
group method (for a detailed description see Ref. l2^ . 
Although the bandwidth is different in the two bands, in 
the logarithmic approximation a common cut-off can be 
used,^^ i.e., in the argument of the logarithm, a quantity 
uj/Eq will be used instead oiuj/Ea, oj/Ei, or 2uj/ (Ea+Eb). 

Up to second order in the couplings the renormaliza- 
tion-group procedure leads to the following scaling equa- 
tions (here x = Eq/Eq, and prime denotes the renormal- 
ized cut-off): for the forward-scattering terms without 
charge transfer between the bands 



dgj_ 

dlnx 



ct\2 



(25a) 
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+ (53^Y)' + (.93T'')' 
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for the backward-scattering terms without charge trans- 
fer 



dgl 
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+ (n - 2) (.g^l)^ + (.gJl)^ + (33^")^ + (53I)' 



for the charge-transfer terms 
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+ 2(n-2).g^l.g3t'^^ + 253tg^l, 
for the intraband umklapp processes 
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and for the interband umklapp processes 
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These systems of equations are the generalizations of 
the set of equations investigated in detail in Ref. to n- 
valued spin and to the case when umklapp processes are 
taken into account. As it will be discussed, depending 
on the filling of the bands some or all umklapp processes 
could be neglected, because momentum conservation can- 
not be satisfied for states near the Fermi points. 

Once the scaling equations have been derived one can 
look for their fixed points. Two types of fixed points 
can be distinguished. One of them corresponds to fixed 
points, where some of the couplings vanish while all oth- 
ers remain weak, the dimensionless couplings introduced 



in H23|l are less than unity. In what follows the term 
fixed point is used even if the marginal couplings scale to 
non-universal values on an extended hypersurface in the 
space of couplings. 

In the other type of fixed points at least some of the 
couplings scale to large values. Using the scaling equa- 
tions derived in the leading logarithmic approximation 
these couplings scale to infinity. When the next-to- 
leading logarithmic corrections are taken into account, 
finite strong-coupling fixed points are obtained, but they 
are located outside the region of validity of the low-order 
calculation. On the basis of this perturbative approach 
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one cannot decide whether these couphngs scale to a fi- 
nite fixed point with dimensionless strength larger than 
unity or to infinity. Fortunately, for the characterization 
of the phases it is sufficient to know whether these cou- 
plings are relevant or irrelevant. 



IV. TOMONAGA-LUTTINGER FIXED POINT 

In the most general case the model is characterized 
by 18 coupling constants, six of them belong to forward- 
scattering processes without charge transfer, three are 
backward-scattering processes without charge transfer, 
three correspond to charge-transfer processes and six to 
umklapp processes. When all true backward-scattering 
processes are irrelevant, their coupling constants scale to 
zero, 



9i± = 9i± 



ab* 
9lJ_ 



0, 



(30) 



all processes in which charge is transferred from one band 
to the other are also irrelevant, 



(31) 



and all umklapp processes are irrelevant, the right-hand 
sides of the scaling equations for the couplings of the 
forward-scattering processes vanish also. This fixed point 
will be denoted as FP~TL, referring to the fact that only 
those processes survive in which charge and spin conser- 
vation are statisfied in each branch of the spectrum, and 
therefore Tomonaga-Luttinger behavior is expected. 



A. Stability analysis 

The behavior of the couplings close to a fixed point 
can be studied analytically using the eigenvalues of the 
linearized matrix 



(32) 



{a*} 



of the renormalization-group transformation, where 
prime denotes the renormalized coupling, i and j are 
short-hand notations for all indices, and * refers to the 
fixed-point value. Deviations from the fixed-point are 
relevant if the eigenvalue is larger than unity, they are ir- 
relevant if the eigenvalue is less than unity. A unit eigen- 
value indicates marginal direction. Whether a marginal 
direction is marginally relevant or marginally irrelevant 
is determined by higher-order corrections. 

First the case of generic band filling is studied, when 
all umklapp processes can be neglected. In the weak- 
coupling regime near FP-TL the matrix Lij defined in 
(|32|l can be diagonalized. The 12 eigenvalues are: 



Ai = l + 2(.gf +52"l)lna;, 
A2=l-l-2(5f|*+.g2^1)lna;, 



(33a) 
(33b) 



A3 = l + 2(.gf* + .g2"i*)lnx, 
A4 = 1 + (-.gf - .gf + 2gf*) \nx, 
A5=l + (.9S+.92l + 2.gf*)lna;, 
A6=l + (.92"l+.92l-2.g2"l*)lna; 



(33c) 
(33d) 
(33e) 
(33f) 



and the six other eigenvalues are equal to 1. The fixed 
point can be attractive if the eigenvalues Ai, . . . , Ae are 
not greater than unity. Since during the renormaliza- 
tion procedure the cut-off is scaled to smaller values, the 
logarithmic multiplicative factor has negative sign, the 
coupling constant combinations appearing in the brack- 
ets have to be non-negative in the irrelevant directions. 
Thus for generic band filling FP-TL can be attractive if 
the following inequalities hold: 



5|| 



> - 



9t> 



52_L' 
.92I, 
.92I* 



9r ^ 

2g^^''*>gr+9- 
52"l+.92l>252t* 



6* \, a* 6* 

^ ^ 52_L 52±> 



(34a) 
(34b) 
(34c) 
(34d) 
(34e) 
(34f) 



Although these relations give the stability condition of 
the fixed point only and cannot be used to get the whole 
basin of attraction, clearly the Tomonaga-Luttinger fixed 
point has an extended basin of attraction in the param- 
eter space. 

When either band a or band b is half filled, i.e., Aka or 
Akh is equal to a reciprocal lattice vector G, the contri- 
butions of the intraband umklapp processes have to be 
taken into account. We will consider the process g^^^ but 
the results can be easily extended to the case when g^j^ is 
relevant by replacing the couplings gf by the correspond- 
ing g^. FP-TL is recovered if besides the requirements 
given in (|3UII (backward-scattering processes are irrele- 
vant) and l|31|l (charge-transfer processes are irrelevant) 
the intraband umklapp process, too, is irrelevant, the 
fixed-point value of its coupling vanishes: 



53I 



0. 



(35) 



With one new coupling added the 13 eigenvalues of 
the stability matrix Lij are as follows: 12 of them are 
the same as given earlier in H33|) . and the 13th eigenvalue 
is 



Ai3 = l + 2(.gr-52°l)lnx. 



(36) 



The condition of stability of the fixed point is similar 
to the case discussed above, except that the requirement 
9«* > 52I changes the first inequality to 



> I.92II 



(37) 



This means that Luttinger-liquid behavior is expected in 
a smaller region of the parameter space. 
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The situation is similar when the system is half-filled, 
2{ka + kb) is equal to a reciprocial lattice vector G and 
therefore the interband umklapp processes have to be 
taken into account. In this case the two Fermi veloci- 
ties are equal (7 = 1). A weak-coupling FP-TL is ob- 
tained if besides the backward-scattering processes 
and charge-transfer processes ()31f) all interband umklapp 
processes are irrelevant, 



ahab^ 
931- 



abba 
931- 



aabb^ 
931- 



am 



0. 



(38) 



In this case 16 coupling constants have to be taken 
into account. 12 of the 16 eigenvalues are the same as for 
generic filling and the remaining 4 are: 



Al3 
Ai4 



1 + 2(.gf * 
1 + (-.92I 



.92T)ln2;, 



.92^1 



2.9y)lna:, 



Ai5=l + (.9r+.9r-2ff2"i*)lna^, 
Ai6=l + (-.gr-.9||*+2.9l[''*)lna;, 



(39a) 
(39b) 
(39c) 
(39d) 



The attractive region of the weak-coupling fixed point is 
determined by the inequalities: 





9r> -92L 


(40a) 




9t> - .92I, 


(40b) 




9t*>\92l*l 


(40c) 




2g||''* >5||* + .gf , 


(40d) 




2<*>|.g2"l+52'll, 


(40e) 


52^1 


+ .92l>2.g2°l*, 


(40f) 


.911* 


+ .9r>2.92"l* 


(40g) 



The basin of attraction of the Tomonaga-Luttinger fixed 
point is further reduced by the requirement that all in- 
terband umklapp processes should be irrelevant. 

One can see that FP-TL cannot be reached when the 
model is fully attractive, when the initial values of all 
couplings are negative. On the other hand in a fully re- 
pulsive model, when the initial values of all couplings are 
positive, Luttinger-liquid behavior is a possibility. Unfor- 
tunately, as will be seen, the repulsive Hubbard ladder 
does not belong to this class. 



Bosonization of the Tomonaga-Luttinger part of 
the model 



It is expected that for such couplings which scale to the 
weak-coupling FP-TL, this n-component fermionic lad- 
der has 2n gapless bosonic excitations, two soft charge 
and 2{n — 1) soft spin modes, and therefore the model 
is equivalent to a 2n-component Luttinger liquid. This 
is shown here by transforming the fixed-point Hamil- 
tonian into bosonic language. For this a more general 
model will be considered where only intraband and in- 
terband forward-scattering terms are allowed without 



charge transfer, but the couplings g2-„„f, 52;ctct'' 

52 CTcr' depend on a and a' separately. It is assumed, 

however, that they are symmetric in the spin indices. 



a{b,ab) a{b^ab) 

92;aa' ~ 92-rj'a 



(41) 



Introducing the density operators corresponding to 
particle-hole excitations with small momentum. 



(9)=i:&L.>d 



(42) 



k 



the kinetic energy of the particles in the two bands given 
in Q can be written as usual in the form 



iThY, / dx\va{nl,{x)Y + v,{nX,{x)Y] , (43) 



where A = ± and the operators ru^'^^ {x) are the densities 

in real space, the Fourier transforms of "■^''ct (?)■ The 
low-energy part of the spectrum contains excitations with 
large momentum, too. They correspond to excitations in 
which the number of particles with a given spin index in 
a given branch changes. They can be taken into account 
by giving the change in the number of particles in the 
branches. For simplicity we will not write out these terms 
since they do not influence the results. 

The interaction terms containing intraband and inter- 
band forward-scattering processes without charge trans- 
fer can also be expressed in terms of the small-momentum 
components of the densities: 



'T/Lutt 

''•int — 



+ Vb92-aa'<Al)^-,'r'i-l) 

+ liva+vt)gt.'nl^Aq)nt^,,{-q) 
which can be written in real-space representation as 



(44) 



n 



Lutt 



crcr' 



da; 



'^aff2;a^'<,a(a;)'^^,^'(a;) 



(45) 



+ «b52;<T<T''^+.<T(2;)?^-,a' {x) 

+ \{va + v,)g^l,,n\^,{x)n'_^,,{x) 

+ \{va + Vb)g^'^^^,n\^^{x)n''^^,,{x) 



First, the Hamiltonian is diagonalized in its spin in- 
dices (7 and a' . To this end we perform an orthogonal 
transformation to a new basis corresponding to such com- 
binations of the bosonic densities, that two of them (one 
for each band) are symmetric and 2(n — 1) {n—1 for each 
band) are antisymmetric in the spin quantum numbers: 



a(b) 



{x) 



1 



x) , 



(46) 
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a{b) 
A,ms 



where m = l,...,n — 1. This transformation loads to 
a Hamiltonian in which the spin and charge degrees of 
freedom are separated, 

n = Y, [ da;(Hki„,,(a;) + Hint j (a;)) , (47) 

i 

where j = c, Is, 2s, . . . , (n — l)s, and the kinetic energy 
part of the Hamiltonian density is 

H^inA^) = TTfi^ [va{nlj{x)Y + v,{nij{x)f] , (48) 



while the interaction term can be written as 



(49) 



Hintj(a;) = 27rh[vag2-jnl.j{x)ntj{x) 

+ «652u"+j(2;)'^-j(a;) 

+ b2-j{va + Vt)nlJx)ntAx) 



+ 



\92ui'^a + Vb)nXj{x)n''_j{x)), 



In this representation, the Hamiltonian can be written 

as 

n = ^ f dx{n'^{x) + n]{x) + nf{x)) , (55) 
j ■' 

where 

n'^.{x) = ^\^^K^,{W,{x)y 
with 



1 



(56) 



1+52-;/ 



(57) 



and similar expressions with index b, and the coupling 
between the bands is given by 



nf{x) 



K>'a + >'h) 



-^^w^{x)n){x) 



(58) 



and the corresponding couplings are 



a(6,afe) 
92;c 



a(b,ab) 
92-,ms 



1 " 

- E 



a{h,ab) 
92;a(T' ' 



(50) 



1 



Ea(b,ab} m a{b,ab) 

92;aa' + ^ _|_ S'2;m+l,m+l 



m(m + l) 

1 / a(b,ab) a{b,ab) \ 

^ ^ X \92;cr,m+l + 92;m+l,a J ■ 



o-=l 



In order to obtain the velocities of the different modes 

the Hamiltonian is rewritten in terms of continuum 
bosonic fields. The fields {x) are related to the den- 
sities by the relation 



imlfix) = -d^(f)l)"/{x). 
The total phase fields and their duals are 

cpf\x) = <i>i^^jix)+r-^^lix), 
9f\x) = r^^{x)-r-^^{x). 

The fields H"'''^ defined by 



Mb), 



(51) 



(52) 



7rnf>{x) = -djf>{a 



(53) 



are canonical conjugates to the fields (j)^^''\x), they sat- 
isfy the commutation relation: 



[Iif\xUf\x')\=^5,,.6{x-x'). 



(54) 



The first two terms of the Hamiltonian can be inter- 
preted as the Hamiltonian of two n-componcnt Luttingcr 
liquids. The coupling 7i"^(x) is diagonal in the compo- 
nent index, it couples the 2n components into n pairs. 
Without this coupling, in the most general case, the 2n 
modes have different velocities, so not only spin-charge 
separation occurs in the system but all modes appear sep- 
arately. In the special case, when the couplings depend 
only on whether a = u' ov a ^ a' , the couplings are 



9^^^=9t^ + {n-l)g<l\ 



(59) 



and 



Mb) _ Mb) _ a{b) 
i/2;ms ~ i/211 i'2_L 



m = l,2,...,(n-l), (60) 



so the n — 1 spin modes in band a (6) have the same 
velocity. 

The full Hamiltonian can be diagonalized by the "sym- 
metric" and "antisymmetric" combinations of the oper- 
ators belonging to the two bands: 



\/Ua,j + Ub,j 



(61) 

and similar expression holds for the conjugated momenta: 



V"a,j + Ub,j 



m ^Ua,jKa,j±m Vub,jKb,j , (62) 



such that they satisfy the canonical commutation rela- 
tion. 
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With these operators the Hamiltonian density is the 
sum of 2n Tomonaga-Luttinger Hamihonians: 



(63) 



where 



(64) 

and the velocities and Luttinger parameters are given by 



K. 



+,0 



(65) 



where 



9j 



52;jK+f&) 



(66) 



The velocities are the same in the two sectors, there- 
fore, in the most general case only n different velocities 
are found for the 2n modes. If the couplings depend 
on the relative spins of the scattered electrons only, 2 
different velocities are found and one recovers the usual 
spin-charge separation. 



V. STRONG-COUPLING FIXED POINTS 

In all other solutions of the scaling equations, some 
of the backward-scattering, charge-transfer or umklapp 
processes are relevant and scale to a strong-coupling fixed 
point. Although these processes may drive the marginal 
forward-scattering processes, too, to the strong-coupling 
limit, these fixed points will be classified by giving those 
backward-scattering, charge-transfer and umklapp pro- 
cesses only, which become relevant. 



A. Fixed points 

First the backward-scattering terms are considered. It 
is easy to check that the scaling equations can be sat- 
isfied if the charge-transfer and umklapp processes are 
still irrelevant, they scale to zero, and only some of 
the backward-scattering processes — together with some 
forward-scattering ones — scale to strong coupling. FP- 
BS-a(b) denotes the fixed point where the coupling of 
the intraband backward scattering in band a (or in band 
b) has non- vanishing value. The fixed point FP BS-ab 
corresponds to the case when the interband backward 
scattering is relevant and the intraband backscatterings 
are irrelevant. It is possible that one or both of the in- 
traband and the interband backward-scattering processes 



are relevant at the same time, but as we will see these 
fixed points do not give a new phase. 

Similarly, it can be checked that the scaling equations 
can be satisfied by assuming that in the fixed point, one of 
the charge-transfer forward-scattering processes has non- 
vanishing coupling, while all bac;kward-scattering and 
umklapp processes are irrelevant. Accordingly FP CT- 
FS|| and FP-CT-FS_L denote the fixed points where ^[j^ 

or g2±, respectively are relevant, while the other charge- 
transfer processes arc irrelevant. 

While the backward-scattering and charge-transfer 
forward-scattering processes may become relevant on 
their own, this is not the case for the charge-transfer 
backward-scattering process gj*,^, except if n = 2. For 
n > 2, if the charge-transfer backward-scattering pro- 
cess is relevant, it will necessarily drive the interband 
backscattering processes {gij_) and the charge-transfer 
processes between electrons of parallel spin isf"), too, to 
strong coupling. The corresponding fixed point is de- 
noted by FP-CT-BS. 

Similar situation occurs when one of the bands is half- 
filled. In this case the intraband umklapp processes may 
become relevant, and, if n > 2 these processes drive the 
intraband backward-scattering processes, too, to strong 
coupling. This fixed point is denoted by FP-U-a(b). 

The names of the fixed points FP-U-abab, FP-U-abba, 
FP-U-ab||, and FP-CT-U, which may be reached when 
the system is half filled, are self-explanatory. E.g., in 
the fixed point FP-U-abab, the coupling ^g'j^"'' scales 
to strong coupling, while in the fixed point FP-U- 
abba the coupling .g"^'"' is relevant. For n = 2, these 
fixed points can be reached without driving charge- 
transfer, backward-scattering or other umklapp processes 
to strong coupling. For n > 2, however, necessar- 
ily some other couplings are also relevant. E.g., the 
coupling ^3^"'' will drive the couplings gj*! ^'^d Ssy to 
strong coupling, or together with g^^'^ both the intra- 
band backscattering and the umklapp process g^ become 
relevant. The charge-transfer umklapp process drives the 
interband backward-scattering processes to the strong- 
coupling limit. 

Depending on the band filling other, more complicated 
fixed points may also be reached if appropriate initial 
couplings are chosen. They correspond to regions where 
several of the backward, charge-transfer or umklapp pro- 
cesses become relevant at the same time, although they 
are not necessarily driven to strong coupling by the other 
relevant processes. 



B. The role of the individual terms 

As we have seen the two-leg fermion ladder with n- 
valued spin has two charge modes and 2(n — 1) spin 
modes. In the strong-coupling fixed points, where one 
or more backward-scattering, charge-transfer or umklapp 
processes are relevant, gap may be opened in the spec- 
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trum of some of the modes. The bosonization proce- 
dure allows us to determine how many of the modes are 
gapped and how many remain gapless. For this we re- 
express the one-particle fermion fields and the relevant 
scattering processes in terms of the boson degrees of free- 
dom. 

The continuum fermion fields of the right and left mov- 
ing particles are defined by 



(67) 



and by similar formulas for the operators b which annihi- 
late left or right moving particles in band b. The bosonic 
phase fields are related to the fermionic fields by the re- 
lations 



(68) 



2tt 



Here F^^'^ are the Klein factors which make sure that the 
fermion operators satisfy the proper anticommutation re- 
lations even though the phase fields 0^*^^ satisfy bosonic 
commutation relations and they commute with the Klein 
factors. The total bosonic fields and their duals are: 



a(b) 

— ,(T 1 



(69) 



These fields are related to the charge and spin fields de- 
fined earlier by the relations 



1 " 



(b) 



Mb) 




(70) 



:0.{b) _ 



Mb) 



with m = 1, . . . , n — 1. 

In order to see what the role of those scattering pro- 
cesses is which should vanish in the Tomonaga-Luttinger 
fixed point and which may drive the system into a strong- 
coupling fixed point, we rewrite them in this continuum 
boson representation. The true backward-scattering pro- 
cesses are proportional to 



9\ 



ah 
9l± 



dx , 



dx , 



Ecos(v^(<^: 



C' - (it')) 



xcos(^^F(e«-C-^'+^'')) 



da; . 



It follows from the sine-Gordon form of the interactions 
that they may give rise to gapped soliton solutions either 
in the charge or spin sector. It is easy to sec that the 
backward-scattering terms contain spin modes only. For 
this one has to realize that the combinations — 0^1''^ 
(and the similar combinations of the dual fields) occur- 
ring in the cosine functions are orthogonal to the charge 



modes 



Mb) (Mb) 



). The intraband backward-scattering 
process opens n — 1 spin gaps while the interband process 
{gfl) opens 2(n — 1) spin gaps. 

In boson representation, except for a factor, the charge- 
transfer processes have the form 



9\\ I Y.^os{s/^{ei-el))\dx, 



xcos(v^(0^-e-</'^ + 



dx , 



(72) 



xcos(v/^(C-€'+<^a-<^^')) 



da; . 



The first two expressions which describe charge-transfer 
forward-scattering processes, contain one charge mode 
and n — 1 spin modes, moreover, they couple the same 
modes, and therefore they modify the spectrum in the 
same way. For this reason we will not distinguish be- 
tween these processes and the notation FP-CT-FS will 
be used for any of the two strong-coupling fixed points 
FS CT-FSIl and FS-CT-FS_L or for the fixed hypersur- 
face FP-[CT-FS||]-F[CT-FS_L] on which both ^^f' and g^]_ 
are relevant. In the expression for the charge-transfer 
backward-processes the first cosine function contains one 
charge mode and n — 1 spin modes if n > 2, while if 
n = 2 it contains only a charge mode. The second cosine 
function contains the other n—1 spin modes. As a result 
one charge mode and one spin mode appears in the 
term if n = 2, but one charge and 2(n — 1) spin modes 
for n > 2. 

The umklapp terms may also be responsible for the 
opening of a gap in the excitation spectrum. The intra- 
band umklapp processes have the form: 



93± 



Ecos(V4;^(C + C')) 



(71) gl^J\j2cos{V4^{cf>l + ^l)) 



dx , 
da; . 



(73) 



These umklapp processes do not contain spin modes but 

one charge modes only if n = 2, while if n > 2 the charge 
mode is coupled to n — 1 spin modes. 
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The interband umklapp processes take the form 
xcos(V^(C + e+0' + 0'')) da:, 

93± J E (^('^- + + + ) 



(74) 



xcos(v^(^?^+e-e^-e^')) 



The interband umklapp scattering g^ involves one 
charge and n — 1 spin modes for any n. For the other 
terms different behavior is found for n — 2 and n > 2. 
The processes 53^"'' and g^^°' involve the same modes, 
namely they open one charge gap and 2{n— 1) spin gaps 
for n > 2, however, only one charge gap and one spin gap 
are opened if n = 2. Therefore similarly to the charge- 
transfer forward scatterings they (alone or together) pro- 
duce the same gap structure in the spectrum. In what 
follows these processes will not be distinguished and the 
notation FP-U-ab_L will be used for any of the fixed 
points FP-U-abab and FP-U-abba or for the hypersur- 
face on which g^'^'' and g^f"^ are relevant simultaneously 
(FP-[U-abab] + [U-abba] ) . 

The term corresponding to charge-transfer umklapp 
processes contains only two charge modes if 71 = 2, but 
two charge and 2{n — 1) spin modes when n > 2, i.e., all 
modes are gapped. 



C. Possible phases 

Following Balents and Fisher^^ the phases belonging 
to the various fixed points will be characterized by the 
number of gapless charge and spin modes, respectively. 

In the attractive regime of the Tomonaga-Luttinger 
fixed point where all backward-scattering, charge- 
transfer, and umklapp processes are irrelevant, their cou- 
plings scale to zero, both charge modes and all spin 
modes are gapless. In this C2S2(n — 1) phase the sys- 
tem is a 2n-component Luttinger liquid. 

As it has been pointed out earlier, the backward- 
scattering processes are responsible for opening gaps in 
the spin sector, while the charge-transfer and umklapp 
processes may open gaps both in the charge and spin 
sectors. Because of the special behavior for n = 2 this 
case will be discussed separately. 

For generic filling and for n > 2 the charge sector re- 
mains gapless if the charge-transfer processes are irrel- 
evant. Depending on whether one of the intraband or 



the interband backward-scattering processes is relevant, 
and the fixed point FS-BS-a(b) or FS-BS-ab is reached, 
a phase C2S(rt — 1) or C2S0 is found, as shown in Ta- 
ble m Since all spin modes are gapped if the interband 
backward-scattering is relevant, the gap structure is not 
modified by the intraband backward-scattering processes. 
The same phase is obtained irrespective whether these 
processes are relevant or not. Similarly, because the two 
intraband backward-scattering terms open gaps in dif- 
ferent spin modes, the spin sector is fully gapped when 
both of them are relevant. The phase is not modified by 
a relevant or irrelevant interband backward scattering. 



fixed point 


pliase 


FP-TL 


C2S2(n - 1) 


FP-BS-a(b) 


C2S(n- 1) 


FP-[BS-a] + [BS-b] 
FP-BS-ab 


C2S0 


FP-CT-FS 


C1S(71 - 1) 


FP~CT-BS 
FP-[CT-FS]-H[BS-a(b,ab)] 


CISO 



TABLE I: Possible fixed points and phases of the n- 
component fermion ladder for n > 2 at generic filling. 

The phase realised in the regime where the charge- 
transfer forward-scattering processes are relevant, can be 
described as ClS(n — 1). Finally if both the backward- 
scattering and charge-transfer processes are relevant, the 
spin sector is fully gapped and only one charge mode is 
gapless (CISO). Due to the fully gapped spin sector in 
the fixed point FP~CT-BS {gl\^ is relevant), the same 
phase is obtained for n > 2 even if the other backward- 
scattering processes become relevant. 

As we see there is always at least one gapless charge 
mode, and a gapped charge mode necessarily implies that 
some of the spin modes become gapped. 

The situation is similar for n = 2, except that if 
the charge-transfer backward-scattering processes are rel- 
evant, they lead to phase ClSl. However, if either 
the charge-transfer backward and forward-scattering pro- 
cesses are relevant simultaneously (in this case the cou- 
plings related to other backward-scattering processes 
may be relevant or irrelevant) or the charge-transfer 
backward-scatterings and at least one of the other 
backward-scatterings are relevant simultaneously, the 
whole spin sector is gapped. The phases corresponding 
to the various fixed points are given in Table ITTl 

When one of the subbands is half filled, and the corre- 
sponding intraband umklapp processes are relevant, the 
coupling (73_L (or g^^) opens one charge gap and n — 1 
spin gaps. If the backward-scattering processes become 
relevant simultaneously, they open gaps in the remaining 
gapless spin modes, while charge-transfer together with 
the umklapp processes opens gaps in all branches of the 
spectrum. This, of course, is not modified if besides the 
umklapp and charge-transfer processes other couplings, 
e.g., backward scatterings become relevant. The possible 
phases are given in Table Ulll 
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fixed point 


plicLse 


FP-TL 


C2S2 


FP-BS-a(b) 


C2S1 


FP-[BS-a]+[BS-b] 
FP-BS-ab 


C2S0 


FP-CT-BS 
FP-CT-FS 


ClSl 


FP-[CT-BS]+[BS-a(b,ab)] 

FP-[CT-BS]+[CT-FS] 
FP-[CT-FS] + [BS-a(b,ab)] 


CISO 



TABLE II: Possible fixed points and phases of the n- 
component fermion ladder for n — 2 sd, generic filling. 



fixed point 


phase 


FP-TL 


C2S2(n - 1) 


FP-U-a 


ClS(n- 1) 


FP-[U-a]+[BS-b(ab)] 


CISO 


FP-[U-a]+[CT-BS] 


coso 


FP-[U-a]+[CT-FS] 



TABLE III: Possible fixed points and phases of the n- 
component fermion ladder with half-filled subband for n > 2. 

The phases obtained for n = 2 are given in Table Hvl 
The difference compared to n > 2 is due to the fact 
that the intraband umklapp processes open gap in the 
charge sector only and therefore phases with a fully gap- 
less spin sector and one soft charge mode (C1S2) and 
a fully gapped charge sector with one soft spin mode 
(COSl) are also possible. Similarly to the n > 2 case 
only one charge gap can be found in the spectrum if the 
charge-transfer processes are irrelevant (phases ClSl and 
CISO). 



fixed point 


phase 


FP-TL 


C2S2 


FP-U-a 


C1S2 


FP-[U-a] + [BS-a(b)] 


ClSl 


FP-[U-a]-h[BS-a]+[BS-b] 
FP-[U-a]-h[BS-ab] 


CISO 


FP-[U-a]+[CT-BS] 
FP-[U-a]+[CT-FS] 


COSl 


FP-[U-a]+[CT-BS]-l-[BS-a(b,ab)] 

FP-[U-a] + [CT-BS] + [CT-FS] 
FP-[U-a] + [CT-FS]-h[BS-a(b,ab)] 


COSO 



TABLE IV: Possible fixed points and phases of the n- 
component fermion ladder with half-filled subband for n = 2. 

The phases occurring when the system is half-filled and 
one or more interband umklapp processes are relevant are 
listed in Table E| for n > 2. The umklapp scattering 
processes — similarly to the charge-transfer processes — 
open at least one charge and n — 1 spin gaps. A soft 
charge mode may be present in the spectrum only if 
all charge-transfer processes are irrelevant, moreover, in 
phase ClS(n — 1) the umklapp scattering between elec- 
trons with parallel spins only may be relevant. As it was 



mentioned in the previous section the couplings g^j^ ^^'^ 
gabba Qpgjj 2(n — 1) Spin gaps, therefore — due to the fully 
gapped spin sector — the relevance of the processes which 
open gap in the spin sector only 5i_li does not 
modify the spectrum. However, in phase CISO, besides 
the coupling at least one backward scattering has to 
be scaled to strong coupling. We note that the charge- 
transfer umklapp process g^^ alone makes the spectrum 
fully gapped, but of course in phase COSO, all couplings 
may be relevant. 



fixed point 


phase 


FP-TL 


C2S2(n - 1) 


FP-U-ab|| 


ClS(n- 1) 


FP-U-ab± 
FP-[U-ab||]-f[BS-a(b,ab)] 


CISO 


FP-CT-U 




FP-[U-ab_L]-h[CT-FS] 
FP-[U-ab_L]V [CT-BS] 
FP-[U-ab||]-(- [CT-FS] 


COSO 


FP-[U-ab||]-|- [CT-BS] 





TABLE V: Possible fixed points and phases of the half-filled 
n-component fermion ladder for n > 2. 

If n = 2 the charge-transfer umklapp processes do not 
open gap in the spin sector, therefore, new phases — 
namely C0S2 and COSl — occur in the phase diagram 
compared to n > 2 (see Table EJ- It has to be noted 
that similarly to n > 2 the spin sector can be fully gapped 
without relevant backward-scattering processes (FP-[U- 
ab_L]-|-[U-ab||]) and the relevance or irrelevance of the 
non-charge-transfer backward scattering processes does 
not modify the gap structure. Moreover, in the fixed 
points, where relevant charge-transfer umklapp processes 
are present, the relevance of any other coupling leads to 
phase COSl (the last three fixed points for phase COSl), 
except for gf]^ which opens 2{n ~ 1) spin gaps, so the 
spectrum becomes fully gapped and the corresponding 
phase is COSO. 



fixed point 


phase 


FP-TL 


C2S2 


FP-U-ab±(ab||) 


ClSl 


FP-[U-ab±]-h[BS-a(b,ab)] 

FP-[U-ab_L]-F[U-ab||] 
FP-[U-ab||]-h[BS-a(b,ab)] 


CISO 


FP-CT-U 


C0S2 


FP-[U-ab±]+[CT-FS] 
FP-[U-ab||]-h[CT-BS] 
FP-[CT-U]-h[BS-a(b,CT)] 
FP-[CT-U] [U-ab± (ab 1 [ )] 
FP-[CT-U]+ [CT-FS] 


COSl 


e.g. FP-[U-ab_L]-h [CT-BS] 


COSO 



TABLE VI: Possible fixed points and phases of the half-filled 
n-component fermion ladder for n — 2. 

It can be seen that all possible phases have been found 
if n = 2 while for n > 2 some theoretically possible phases 
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do not occur. This is due to the fact mentioned above 
that in this case there are no couphngs which open gap 
in the charge sector only: a charge gap opens always to- 
gether with spin gaps. A more complete phase diagram 
could be obtained if the couplings have a more general 
spin dependence. Although, even in this case the charge 
gaps are coupled to spin gaps, it would be possible to 
open a single spin gap — if the corresponding fixed point 
exists — contrary to case investigated above where n — 1 
spin gaps occur always together. 

VI. NUMERICAL RESULTS 

The behavior close to the weak-coupling fixed point has 
been analyzed in Sec. IV analytically. A similar calcula- 
tion close to the strong-coupling fixed points is not pos- 
sible since — as mentioned earlier — the scaling equations 
derived in the leading logarithmic approximation are not 
valid in the strong-coupling regime. The relevance or 
irrelevance of the couplings could, however, be deduced 
from the numerical solution of the scaling equations. We 
have done numerical studies for different band fillings 
for models in which the bare couplings correspond to a 
repulsive Hubbard ladder, and also for somewhat more 
general situations characterized by several bare coupling 
constants. 




5 10 5 10 



-ln[Q/Eo] -Inlto/Eo] 

n=2 n=4 

n=3 11=6 

FIG. 6: The second-order scaling curves for a weak-coupling 
Hubbard ladder {U/{2nhVa) = 0.1, vt/va = 1.5) for the 8 
independent couplings and for different values of n. 



A. Hubbard ladder 

In a Hubbard ladder the initial couplings of all intra- 
band and interband scattering processes gff^^ are related 
to the same on-site interaction U defined in Eq. while 
^ab g^^,^ equal to zcro. The numerical results 

presented below were obtained for U/{2nhva) = 0.1 and 
different values for the ratio Vb/va- Other choices for 
U / {2TThva) gave qualitatively identical results. 

1. Generic filling 

When the umklapp processes may be neglected 12 scat- 
tering processes have to be taken into account. Looking 
at the lowest-order scaling equations it can be shown that 
if the initial couplings statisfy the relations 

a{b.ab, ct) a{b.ab, ct) a{b.ab, ct) 

what is the case for the Hubbard ladder, the same rela- 
tions remain valid for the renormalized couplings. Due to 
these scale-invariant surfaces eight independent couplings 
describe the system. The second-order scaling curves for 
these couplings for v^/va — 1.5 are shown in Fig. El 

One can see that all scaling trajectories go to infin- 
ity so all couplings are relevant for all value of n. The 
weak-coupling fixed point, which is stable for the n = 2 
Hubbard chain, is not reached for the Hubbard ladder. 



Since the backward-scattering and charge-transfer pro- 
cesses are relevant the corresponding phase is CISO. For 
generic filling the SU(ri) Hubbard ladder is not a Lut- 
tinger liquid. 

Similar results are obtained for other 7 values pro- 
vided they are of the order of unity. Solving the scal- 
ing equations formally for larger values of 7, it was 
found that successively some of the charge-transfer and 
backward-scattering processes become irrelevant and new 
phases occur. Namely, if v^/va or Va/vh (since 7 takes 
the same value when Vb/va is replaced by Va/vb) is ap- 
proximately 8.3, all charge-transfer processes, the inter- 
band backward scattering and one of the intraband 
backward-scattering terms g"_L 9i±^ belonging to the 
band with smaller velocity, become irrelevant and FP- 
BS-b or FP-BS-a is reached. The corresponding phase is 
C2S{n— 1). The surviving intraband backscattering 
or g^j^ becomes irrelevant above an even higher value — 
approximately 8.7 — of the velocity ratio Vb/va or Va/vb- 
None of the couplings that could open gap in the spec- 
trum seem to be relevant and the Luttinger liquid phase 
C2S2(rt — 1) would be obtained. These boundaries vary 
little with n. These phases appear, however, in such 
situations where one of the subbands is almost empty 
or almost full, the curvature of the dispersion relation 
may be important, and a model with linearized spectrum 
may not be a good approximation. For this reason these 
phases are not indicated in the phase diagram. 



15 



S. Role of intrahand umklapp processes 



If Aha is equal to a reciprocal lattice vector G and the 
umklapp process g^j_ is relevant the surfaces determined 
by H75() remain invariant. Therefore out of the 13 cou- 
plings which give contribution 9 are independent. The 
scaling curves for these couplings obtained in the leading 
logarithmic approximation can be seen in Fig. \7\ 




^ I n=4 

2 4 n=6 

-Inlto/EflJ 



FIG. 7: The second-order scaling curves for the Hubbard lad- 
der {U/{2-KhVa) ~ 0.1, Vb/va ~ 1.5) with half-filled subband. 



It is found that starting from the bare couplings char- 
acterizing the Hubbard ladder the weak-coupling fixed 
point cannot be reached for any value of n and for 
any value of 7. The strong-coupling fixed point FP-[U- 
a]-|-[CT-FS] where the umklapp processes are relevant 
together with the charge-transfer ones but all backward- 
scattering processes are irrelevant seems to be accessible 
for the Hubbard ladder if n = 2 (independently of the 
value of 7). The corresponding phase is COSl. As has 
been seen above, this fixed point could not be reached 
for n > 2 and phase COSO is realized there. Similar re- 
sult has been found^ for the ID SU(n) chain where it 
was shown that the umklapp scattering suppresses the 
normal backward-scattering processes for n = 2 only. 



3. Half-filled system: interband umklapp processes 

When the whole system is half-filled, the contribution 
of the interband umklapp processes should be taken into 
account. The scaling equations are invariant under the 
interchange of bands a and b 

at = .9? (76) 

and therefore out of the 16 coupling which give contri- 
bution in a half- filled system, 13 are independent. The 
scaling curves for these couplings are shown in Fig. |S| 




n=6 

-ln[(M/Eo] 



FIG. 8: The second-order scaling curves for the half-filled 
Hubbard ladder {U/{2-Khva) = 0.1). 

The weak-coupling fixed point determined by H30|l . 
(|31|l and H38|) cannot be reached starting from the Hub- 
bard ladder. Similarly to the case of generic filling, for 
any n and 7 all couplings go to infinity. But since now the 
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umklapp processes, too, are relevant, the corresponding 
phase is COSO. 

Figure El shows the phase diagram of the Hubbard 
modeL The control parameter is the filling nijn where 
Ui is the particle number per site, nijn = and 1 de- 
note the totally empty and full system, respectively, 6o 
denotes the filling when the upper band is almost empty, 
similarly a\ denotes the filling when the lower band is 
almost full. Close to these fillings the linearization of 
the spectrum is not valid, and the shaded regime around 
these points indicates that the results presented in this 
paper using a linearized spectrum may not be reliable at 
these fillings. For this reason the phases found in these 
regimes are not indicated in the figure. 

COSO COSO COSO 
^ ClS(w-l) I C1S0 | CISO I CISQ I CISO | ClS(ra-l) ^ 



2-L 







C0S1 



1/2 
a) 



COSO 



njn 



COSl 



ClSl I CISO I CISO I CISO I CISO I' ClSl [ 



nJ2 



1/2 
b) 



FIG. 9: The phase diagram of the Hubbard ladder: a) for 
n > 2 and b) for n — 2. 

The upper figure (a) shows the phase diagram of the n- 
component Hubbard ladder for n > 2. Figure (b) is valid 
for 2-component fermions. As one can see the n = 2 
and n > 2 models behave differently in the presence of 
intraband umklapp processes. 



B. Chains with forward scattering only 

As another special case we consider two chains with 
nonchiral intrachain forward scattering only. Although 
the interband scattering processes that arise when the 
one-particle transverse hopping is switched on are all of 
forward-scattering type, the charge-transfer processes do 
not conserve the chiral charge and spin. For this reason 
in general, such a ladder does not behave like a Luttinger 
liquid. 

The couplings which play a role in this coupled sys- 
tem are related to scattering processes between electrons 
of the same spin: giy , gy, g^^ and gy*, and to forward 
scattering processes between electrons of different spins: 
52I ^^"i 92±- The scaling equations for the 
intraband scatterings are: 



din a; 



7(^/r)^ 



(77a) 



din a; 



= 7(32^) 



(77b) 



Similar equations hold for the couplings in band b. 
The Lie equations for the interband scatterings without 
charge transfer are: 



d.9||'' 
din a; 

Mi 

dlnx 



(^/2i)^ 



and for the charge-transfer processes 



d.gf 

dlnx 
dgP' 



= - 9\\ 



.9 - 25y 



dlnx 



52\+52\- 252l 



(78a) 
(78b) 

(79a) 
(79b) 



As can be seen, the coupling constants of the scattering 
processes between electrons of the same spin and of dif- 
ferent spins are not coupled, the two scattering channels 
are separated, and analogous equations are obtained. 

These equations have one weak-coupling fixed point 
only, namely when the processes which transfer charge 
from one band to the other are irrelevant. This fixed 
"point" is a six-dimensional hypersurface in the parame- 
ter space and it can be attractive if 



a* 
52_L 



Sgf * > .gf -\ 

+ 52_L > 2g2_L 



fc* 

9\\ 



(80a) 
(80b) 



As we have seen in Sec. II VI at this fixed point the system 
is a Luttinger liquid. One finds, however, different behav- 
ior depending on whether the interband scatterings are 
scaled out of the problem or not. If these processes van- 
ish, the two charge modes become independent, and simi- 
larly the spin modes are decoupled, while if the interchain 
couplings survive, the charge modes become coupled and 
have identical velocities. Similar situation occurs for the 
spin modes as well. 

The scaling equations have been solved numerically for 
one-parameter models, i.e., for models characterized by 
a single bare coupling constant U . According to (|23|l the 
dimensionless couplings g^f^"*^ appearing in the scaling 
equations are U = 2U /TTh{va + v/3 + + vg). 

We have found, starting from a one-parameter 
Hubbard-like model — the Coulomb interaction between 
electrons of the same spin is zero — that the charge- 
transfer processes become relevant, the stable fixed point 
is FP-CT-FS both for fully repulsive and fully attractive 
models. In fact the sign of the couplings gjj* and 52*1 
not relevant in this respect, so this strong-coupling fixed 
point is reached if 



and 



9 



a(0) 
92± 



a(0) 



6(0) 



a6(0) 



6(0) 
92± 



afc(O) 
92± 



±9T 



0, 



±U. 



(81) 



(82) 



17 



One charge and n — 1 spin gaps open in the spectrum, 
the phase is ClS(n-l). 

Due to the analogous structure of the Lie-equations 
hsted above the same behavior is found if the bare cou- 
phngs g'^^^ are all repulsive or attractive while the bare 

g^2± couplings vanish: 

o(0) 6(0) afc(O) , ct(0) ,fj ,oo^ 

9\\ = 9\\ = .9|| = ±.9|| ' ' = ±U (83) 

and 

q(0) _ 6(0) _ ab{Q) _ , ct(0) _ „ ,x 
92± - 92± - 92± - ±32_L - 

Moreover, if the couplings ^ii"' or g^*^ are all repulsive or 
attractive, the gap structure is not modified when some 
of the other bare couplings appear with opposite sign. 

Luttinger liquid behavior could, however, be found 
when the intraband forward scatterings {g2± and 52 j_) 
are repulsive and the interband forward scattering {g2±) 
is attractive, 

5i°^=0, gT=9T-~9T-±9T-U, (85) 

since in this case the charge-transfer processes scale to 
zero. As we have seen in Sec. IV two different velocities 
appear in the system. Similar situation is encountered 
if the intraband forward scattering between particles of 
the same spin is attractive and the intraband forward 
scattering is repulsive, 

„(0) _ r, aiO) _ 6(0) _ a6(0) _ , ct(0) _ 

9± -0,5211 ~ 92\\ ~ ^92\\ - ±.92_L --l^- (00) 

These results are independent of the value of n. Al- 
though the Luttinger-liquid behavior with four different 
velocities is theoretically possible, for the parameter val- 
ues chosen no such behavior has been found, because 
even though the charge-transfer processes are scaled out, 
the interband forward-scattering processes ^j^^ and g^l 
remain finite. 



VII. CONCLUSIONS 

In the present paper, we have investigated a fermion 
ladder in which two n-component fermion chains are 
coupled by interchain hopping. Using the multiplica- 
tive renormalization-group procedure we have derived 
the second order scaling equations, determined the fixed 
points and analyzed the possible phases. 

One weak-coupling fixed point has been found where 
the system behaves like a Luttinger liquid. In this fixed 
point, the processes which violate the conservation of 
the particle number and spin on each branch separately, 
namely the charge-transfer processes from one band to 
the other, the backward and the umklapp processes are 
irrelevant. In the basin of attraction of the Tomonaga- 
Luttinger fixed point, where all modes, namely the two 
charge and 2(n — 1) spin modes are soft, the Hamiltonian 



has been diagonalized for general spin-dependence of the 
couplings. In this fixed point, different behavior has been 
found depending on the relevance or irrelevance of the 
interband forward scatterings. If they are irrelevant the 
system behaves like two independent Luttinger liquid, 2n 
different velocities characterize the ladder, meaning full 
mode separation. On the other hand, if the interband 
scatterings are relevant, this coupling allows for n differ- 
ent velocities only. In the special case when the coupling 
strength depends on the relative orientation of the spins 
of the scattered particles, the velocities of all spin modes 
become equal, and the usual charge-spin separation is 
recovered. 

Besides this weak-coupling fixed point there are several 
strong-coupling fixed points. The relationship between 
them and the appropriate phases have been determined 
by the bosonization treatment. The possible phases for 
generic filling as well as for half-filled subbands and for 
a half-filled ladder have been described. Since the exact 
location of the strong-coupling fixed points cannot be de- 
termined in the approximation used, the strong-coupling 
regime has been investigated numerically. The scaling 
equations have been solved for some special initial values 
of the couplings, among others for a Hubbard ladder. 

We have found that the generically filled Hubbard lad- 
der has fully gapped spin sector and one gapped charge 
mode (CISO), the system does not show Luttinger liquid 
behavior for any n. The half-filled Hubbard ladder is not 
Luttinger liquid, either, it has not only fully gapped spin 
sector but the charge sector, too, is gapped. The phase 
is COSO. These results do not depend on the initial value 
of the Hubbard U, any nonzero value of the Coulomb 
interaction leads to the same phases. 

The 71 = 2 case, the SU(2) Hubbard ladder has been 
investigated in detailed using similar methods in Refs. 
il^ and .13i . The results agree in that regime where the 
linearized spectrum is a reasonable approximation. Six 
of the seven phases predicted in these papers were found 
by us. The remaining one — as well as two of these six 
phases — are in that parameter range where our approx- 
imation is questionable. Interestingly five of these six 
phases occur in the phase diagram of the SU(n) Hub- 
bard ladder, too. 

The Hubbard ladder with half-filled subband shows 
special behavior in the case n = 2. We have found that in 
this case the backward scatterings are irrelevant, there- 
fore, there is one soft spin mode in the system (COSl). 
On the other hand if n > 2 the situation is the same as 
for a half-filled ladder: all couplings are relevant and the 
corresponding phase is COSO, i.e. all modes are gapped. 
Similar result has been found in our earlier work, for 
the SU(rt) symmetric half-filled Hubbard chain: the half- 
filled system has gapped charge mode and soft spin mode 
(phase COSl) for n — 2, while for n > 2 the system is fully 
gapped (phase COSO). In this sense similar behavior has 
been found for the weakly coupled Hubbard ladder with 
half-filled subband and for the half-filled Hubbard chain. 
This result seems to support the claim^ that the intra- 
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band umklapp processes lead to 'confinement', i.e., to the 
vanishing of the interchain hopping in a two-component 
weakly coupled two-leg Hubbard ladder. 

Finally we have shown that if two n-componcnt Lut- 
tinger chains with repulsive forward scattering are cou- 
pled by interchain hopping then — similarly to the case 
of the coupled generically filled n-component Hubbard 
chains — the ID Luttinger liquid state is unstable not only 
for usual electrons with spin-1/2 but for n-component 
fermions, too. Luttinger liquid behavior could be ob- 
tained if some but not all of the couplings are attractive. 



In the present paper, the gapped phases have been 
characterized by the number of gapped modes only. The 
study of the correlation functions could reveal further 
details of the obtained phases. 
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